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Abstract 

In general relativistic framework of large scale structure formation theory in 
the universe, we investigate relation between density contrast and peculiar veloc- 
ity in the Newtonian gauge. According to the gauge-invariant property of the 
energy-momentum tensor in the this gauge, the velocity perturbation behaves as 
the Newtonian peculiar velocity. In this framework, relation between peculiar veloc- 
ity and density contrast with respect to the Newtonian Peebles formula has an extra 
correction term which is ignorable for the small scales structures. The relativistic 
correction of peculiar velocity for the structures with the extension of hundred mega 
parsec is about few percent which is smaller than the accuracy of the recent peculiar 
velocity measurements. We also study CMB anisotropies due to the Doppler effect 
in the Newtonian gauge comparing with using the Newtonian gravity. 



1 Introduction 

In the large structure formation theory, structures originate from small density fluctua- 
tions that are amplified by gravitational instability. These initial fluctuations are assumed 
to be generated during inflationary epoch and are inflated to the beyond of the horizon 
while the size of horizon remains constant during the inflation. These quantum fluctua- 
tions beyond the horizon freeze and evolve like classical perturbations. They reenter the 
horizon at a later time when the horizon grows to the size of perturbation. The observa- 
tional evidence for these small density fluctuations can be seen at the decoupling epoch 
as the anisotropy of CMB and also existence of the large scale structures in the universe. 
One of the consequences of density fluctuation is metric perturbations in the Friedman- 
Robert son- Walker (FRW) universe and deviation of velocity field from the Hubble law 
which is called peculiar velocity. Hence studying peculiar velocity could be a useful indirect 
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method to find out the structures in the FRW universe. In the Newtonian hnear struc- 
ture formation theory, relation between density contrast and peculiar velocity has been 
obtained by Peebles [|l|. From the experimental point of view, Bertschinger and Dekel 
introduced POTENT method to reconstruct velocity potential from the radial component 
of velocity field . Recently Branchini et al combined measurement of radial velocity of 
galaxies and independent measurement of density contrast to estimate |]^. 
On the scales well below than the Hubble radius, the Newtonian theory of gravitation is 
a good approximation and Peebles' approach is widely used. However for observations 
and simulations are probing scales which are a significant fraction of the Hubble radius, 
the light-cone effect should taken into account Our aim in this work is to generalize 
relation between density contrast and peculiar velocity to general relativistic framework. 
The relativistic linear perturbation of FRW universe can be studied by two type of formu- 
lations exist in the literature. The 'gauge fixing' formulation, originated by Lifshitz 
considers perturbed components of the metric which are related to the energy momentum 
tensor. According to the gauge freedom problem in the theory, it is necessary to put 
a constraint on perturbed part of metric for gauge fixing (e.g. synchronous gauge). In 
the Second approach, initiated by Bardeen in 1980 gauge invariant variables can be 
made by the combining perturbed metric elements. These gauge-invariant variables are 
well known from other physical theories. For example, in classical electrodynamics it is 
usually more physical to work in terms of gauge-invariant electromagnetic fields rather 
than in terms of gauge- dependent scalar and vector potentials. The subsequent works in 
FRW perturbation theory can be found in some textbooks [|7|-110 . 



General relativistic analysis of peculiar velocity has been studied in the 'covariant' fluid 



flow |]TT|, 'Harmonic gauge' |Q and ' quasi- Newtonian gauge' fixing methods. In this work 
we are going to use conformal Newtonian gauge also known as the (Longitudinal gauge) to 
obtain relation between density contrast and peculiar velocity. The conformal Newtonian 
gauge was introduced by Mukhanov et al [0 as a simple gauge used for scalar mode of 
metric perturbation. In this gauge perturbation of metric elements are gauge-invariant 
variables. For the case of absence of non-diagonal space-space component in energy- 
momentum tensor, perturbation of metric can be interpreted as relativistic generalization 
of Newtonian gravitational potential. We generalize this idea for peculiar velocity field 
according to its gauge-invariant property in the Energy-Momentum tensor. Here, we 
restrict ourselves to the fiat universe, in the matter and radiation dominant epochs, and 
obtain relativistic relation between density contrast and peculiar velocity field. It is shown 
that in small scales compare to the size of horizon, where general relativistic effect due 
to the light cone effect is negligible, our formulation reduces to Peebles equation. For the 
structures with 300Mpc extension, the relativistic correction is about one percent. We 
also calculate the relativistic Doppler effect contribution on the CMB anisotropics and 
compare it with the Newtonian one. 

The organization of the paper is as follows. In Sec. II we write down perturbation of met- 
ric in Newtonian gauge. In Sec. HI we introduce energy-momentum tensor and Einstein 
equations in this gauge and In Sec. IV, we obtain relativistic relation between density 
contrast and peculiar velocity for fiat universe in radiation and matter dominant epochs. 
In Sec. V , we obtain signature of Doppler effect on the anisotropics of CMB, using 
relativistic peculiar velocity. We conclude in Sec. VI with a brief summary and some 
discussions. 
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2 Perturbation of FRW in the Newtonian Gauge 



Consider a small perturbation of metric /i^i, with respect to FRW background g^,y 

ds^ = {g^i, + h^^)dx'^dx'^ (1) 

where Qf^i, = air^^dr'^ —jijdx^dx^) and r is conformal time. Greek and Roman letters go 
from to 3 and 1 to 3, respectively. The metric perturbation h^^i, could be categorized into 
three distinct types like Scalar, vector and tensor perturbation. This classification refers 
to the way in which h^^ transforms under three-space coordinate transformation on the 
constant time hyper- surface. A covariant description of tensor decomposition has been 
shown in |13|]. Vector and tensor modes exhibit no instability. Vector perturbation decays 



kinematically in an expanding universe and tensor perturbation leads to gravitational 
waves which do not couple to energy density and pressure inhomogeneities. However, 
scalar perturbation may lead to growing inhomogeneities of matter. 

The most general form of the scalar metric perturbation is constructed by using four 
scalar quantities 0, tp, B and E: 

his) ^ ( 20 -B, \ 



\^ B., 2(^7,, - E.,jj J 

One of the main difficulties of the relativistic structure formation theory is that there is 
gauge freedom in the theory and one can make an artificial perturbation by coordinate 
transformation in such a way that infinitesimal space-time distance between the two event 
remains constant. One way to overcome this problem is using gauge fixing and the other 
way is gauge invariant method 0. In the latter method, some quantities are found to 
be gauge invariant under coordinate transformation. It seems that the gauge-invariant 
quantities are similar to the electric and magnetic fields in the theory of electrodynamics 
where they are measurable quantities in contrast to the vector and scalar potentials that 
are gauge dependent parameters. The simplest gauge-invariant quantities from the linear 
combination of 0,^/^,-8 and E which span the two-dimension space of gauge-invariant 
variables are: 



= 0+^[(fi-^')a]' (3) 
= i)--{B-E!) (4) 



a 
a 
a 

where prime represents derivation with respect to conformal time and index (gi) stands for 
gauge invariance. The above variables were first introduced by Bardeen [^]. Newtonian 



gauge is defined by choosing B = E = . It is seen that in this gauge and -0 become 



gauge-invariant variables and the metric can be written as follows: 

ds^ = a(r)^ [(1 + 20)rfr2 - (1 - 2ip)-fijdx'dx^] . (5) 

3 Einstein Equation in The Newtonian Gauge 

In the perturbed background of the FRW metric, the linear perturbation of Einstein 
equation can be written as follows: 

SG^'^ = STrGST^"^. (6) 
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The left hand side of eq.(^) can be obtained from perturbation of metric in the Newtonian 
gauge and the right hand side of it is obtained by perturbed energy-momentum tensor. 
Restricting our attention to scalar perturbation of energy-momentum tensor, we can 
express the most general first order energy-momentum tensor in terms of four scalars as 
follows 0: 

y [p + p)aui -dp + a-ij J 

where 6p and 6p are the perturbations of energy density and pressure, W is potential 
for the three-velocity field = a(t)u*(a;,r) = V*VV and a is anisotropic stress. Bardeen 
has shown that in the Newtonian gauge: 5T^^^^'> = 5T^y [Q. According to this result, 
= a{t)u^ may be regarded as Newtonian peculiar velocity. For a perfect fluid, the 
anisotropic stress which leads to non-diagonal space-space component of the energy- 
momentum tensor vanishes. In this case it can be shown that (p = ip. Using the metric of 
eq.(ll) in the left hand side and eq.(0) in the right hand side of eq.(^, Einstein equation 
can be rewritten: 

V^^ - 3n(j)' - 3(7^2 - K)(f) = 4(j)G6p, (8) 
{a(j))'. = -4(PG{p + p)u„ (9) 
(j)" + mcj)' + {27{ - K)(f) = 4(j)G6p, (10) 

where Ti is the Hubble constant in conformal time. For simplicity in calculation we 
consider the following change of variable. 

= 47tG{p + pf'^uj = {A-RGf/^H^ -n' + K)a^f^uj. (11) 

By substituting eq.(|ll|) into eq.(^), this equation is rewritten in terms of u as follows: 

6" 

u" - Cs^V'^uj - yu; = 0, (12) 

where 

e = \^f\i-^^Y^'. (13) 

a p + p QTiGpa'^ 
and for the adiabatic perturbations = For the case of matter dominant epoch one 
can obtain exact solution for differential equation (|12]), considering {cg = 0), yields: 

/(It 
w ^^^^ 

Our aim is to obtain density contrast and peculiar velocity in terms of u. Density contrast 
6 = Y can be obtained by dividing eq.(H) to component of the FRW equation 
(7^9 + «:= 8zi^): 

iV^4>-3n^' -3(n^ - k)4>). (15) 



3(7^2 + 

In what follows, Eqs.(^ and ( [T3| ) will be our main equations in the following sections. It 
is seen that these two equations are as function of which, can be obtained by solving 
sq-dUD- In the next section we find an explicit relation between density contrast and 
peculiar velocity in the framework of general relativistic perturbation theory in the New- 
tonian gauge, we restrict ourselves to the case of radiation and matter dominant epochs 
in the flat universe. 
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4 Density Contrast— Peculiar Velocity Relation 



From the recent CMB and Supernova type I experiments, it seems that the curvature of 
the universe is flat (k = 0). In this section according to sequence of the radiation 



and matter dominant epochs from early universe, we apply Eqs. (P) and ([15D for these two 
regimes. 

4.1 Radiation dominant epoch 

In the radiation dominant epoch, for the flat universe, according to FRW equations, scale 
factor evolves like a = a^T. In this case eq.(|l^) reduce to: 

For simplicity in the calculation we consider structures larger than Hubble radius, so the 
second term in the right hand side of eq.([T^) can be ignored with respect to other terms 
and uj{x, r) can be obtained form eq.([l4|). Substituting eq. ([l6D into eq. (|T^ , the dynamics 
of oj is obtained: 

where D{x,t) = sin(z/r) + i?2 cos(i/r)]e*^*'^'^'^ and u = Substituting eq . ([T7|) into 
eqs. (pJ]) , (|1^) and we obtain following equations for r), 6{x,t) and Ui{x,T): 

rhfT T^ = ^ [(z/rcos(z/r) - sin(z/r))Ci - (z/r sin(z/r) + cos(z/r))C6] e^*^'^'\ (18) 



r) = i:^\{\{yrf - 1] sin(z/r) + vr\\ - -{vrf] co^{vt)Cx 

1 



+ ([1 - {vrf\ cos(z/r) + vr\\ - -{yTf\ sin(z/r))C2)]e(*^), (19) 

V-v = ^r— ([(z/r cos(z/r) - cos(z/r))Ci - (z/r sin(z/r) + cos(z/r))C2] e^*''''^^ 

— — — [(z/cos(z/r) — z/^rsin(z/r) — z/ cos(z/r))Ci 

- {usm{uT) + z/Vcos(z/r) - z/ sin(z/r))C2])e(^'=-^\ (20) 

According to our consideration for the structures larger than Hubble radius << 1, 
eqs.(^) and (|20|) are reduced as follows: 

V-^ = -^e(^'=-^)C2, (21) 

S = i-e(*'=-")C2. (22) 



Now one can divide eq . (pi|) by eq . (|22D to obtain the explicit relation between density 
contrast and the peculiar velocity: 

V -v = ~H6{uTy. (23) 

It is seen from eq.(^) that in the radiation dominant epoch, the relativistic peculiar 
velocity for a given density contrast is much smaller than what one can expect from the 
Newtonian formula. 
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4.2 Matter dominant epoch 

In the matter dominant epoch, for k = 0, eq.(|l^) is simphfied into this form: 



e = -. 

a 



(24) 



For the flat universe, the scale factor grows as a = aoT^/2. By substituting eq.(|^) in 
eq.(|l^, the dynamics of as a function of conformal time is obtained: 



(25) 



where £"1(0;) and E2{x) are arbitrary functions of spatial coordinates. Using u in the 
eq.(|Tl]), the value of (p obtain as follows: 



{x,T)=Ciix)+T-'C2ix). 



(26) 



Substituting eq.(PBD into eqs.(^ and (|T3|), 6 and u are obtained as functions of conformal 
time, Ci(x) and C2(x) as: 



UAX,T 



^ {V'C,{x) - 81C2{x)) - T-'{r'V'C2{x) + 18^2(0;)) 
1 ,C2(x), 2Ci(x).i, 



as r" 



3 r 



Neglecting the decaying modes, we rewrite eqs.(P7D and (|28|) in the Fourier space: 

1 



4 

ik ■ Vi{k) 



-l{k'T' + l2)C,{k), 



_2_^ 

3ai T 



Ci{k). 



(27) 
(28) 

(29) 
(30) 



Dividing eq.(^0|) by eq.(^), one can obtain an explicit relation between density contrast 
and peculiar velocity in the following form: 



ik.Vk 



fc2 



(31) 



where, H is the Hubble constant in physical time and Vk is the Fourier transformation of 
the peculiar velocity. By using the definition of scalar potential for the peculiar velocity, 
Vk = ^yVk, eq.(0) can be written: 



H6, 



(32) 



In the real space, eq.(P2|) changes to: 



V-v = -H"- + 3if2>V 
h 



(33) 



where h is biasing factor and W = J v ■ dl is the scalar potential of velocity field and can 
be obtained by radial integrating of velocity field. The second term in the right hand side 
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of eq.(^) is a result of general relativistic corrections. This correction can be estimated 
by dividing the second term by the first term in the right hand side of eq. (|32D . 



where dn is the size of horizon and A = | is the size of structure. Is is seen that for the 
small scale structures ( compare to the size of horizon), the general relativistic correction is 
negligible. This correction for the structures with SOOMpc extension is about one percent. 



5 The Effect of Relativistic Peculiar Velocity on CMB 

The inhomogeneities in the universe induce anisotropics in the distribution of relic back- 
ground of the photons on the CMB. The anisotropy of CMB can be caused by several 
reasons. One of the main reasons is that the matter which scattered the radiation in 
our direction had a peculiar velocity with respect to comoving frame when the scattering 
occurred. Since the peculiar velocity of matter is different at different directions in the 
sky, this leads to an anisotropy on the CMB. The universe at the last scattering epoch has 
a redshift of z ^ 1000. Since the equality of matter and radiation occurred at z ^ 5 x 10^, 
we can consider that at last scattering epoch universe resides in the matter dominant. 
Using eq.(|3lD for relativistic peculiar velocity in matter dominant epoch, the relativistic 
anisotropy due to Doppler effect on CMB is obtained as follows: 



5T ^ 



A 



-5 (35) 



where A is the size of structure and dn is the Hubble radius at the decoupling epoch. It is 



seen from eq. (|35|) that the temperature perturbation of the CMB due to peculiar velocity 



in the relativistic framework can be expressed in the term of Newtonian one: 



^ST, (f) 

1 + 3(^^2 



Rel - "T— 777X^2 

where (^) New = is the contribution of Newtonian Doppler anisotropy, 9 = 3A" {Qh){^^^) 
is the angular size of structure M and 9h — 1° is the angular size of Hubble radius. It is 



seen from eq.(|36|) that for the large angles, the relativistic contribution of Doppler effect 
due to the peculiar velocity on the fluctuations of CMB is smaller than the contribution 
of Newtonian anisotropy. 

We compare the relativistic Doppler effect also with Sachs-Wolf effect (which is one of 
the reasons for the the anisotropy of CMB arising from the variation in the gravitational 
potential at the last scattering [1^). It is shown that the contribution of this effect on 



anisotropy of CMB can be given by: 

-T = 3^^) ^^^^ 

Eqs.(|35|) and (|37| ) shows that Sachs-wolf effect at large angles {9 > 9h) is dominant term 
on the anisotropics of CMB while for the small angles Doppler term is dominant. Fig. 
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1 shows the fluctuation of temperature as a function of angular size of the anisotropy 
for the case of Newtonian and relativistic peculiar velocity and Sachs-wolf effect which 
normalized to the density contrast. It can be shown that the contribution of relativistic 
peculiar velocity affect the power spectrum of CMB. 

6 Conclusion 

In this letter, we have tried to identify the peculiar velocity in general relativistic theory of 
structure formation. We obtain a relation between the density contrast and the peculiar 
velocity for the matter and radiation dominant epochs. It has been shown that the 
relativistic correction is about few percent for structures with extension of hundred mega 
parsec. The effect of the relativistic peculiar velocity as Doppler effect on the CMB also 
has been calculated. It was shown that the contribution of relativistic peculiar velocity 
on the anisotropy of CMB is less than what one expected from the Newtonian theory. 
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7 Figure Caption 

This figure shows the normahzed perturbation of temperature to density contrast {^)/S 
as a function of normahzed angular size of structures to the size of horizon j^. Sohd hue, 
dashed hne and dashed-dot hne represent the contribution of Newtonian, Relativistic 
Doppler effect and Sachs- Wolf effect on the anisotropy of CMB, respectively. 
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